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Abstract 

We show that the dispersive force between a spherical nanoparticle (with a radius < 100 nm) and 
a substrate is enhanced by several orders of magnitude when the sphere is near to the substrate. We 
calculate exactly the dispersive force in the non-retarded limit by incorporating the contributions 
to the interaction from of all the multipolar electromagnetic modes. We show that as the sphere 
approaches the substrate, the fluctuations of the electromagnetic field, induced by the vacuum and 
the presence of the substrate, the dispersive force is enhanced by orders of magnitude. We discuss 
this effect as a function of the size of the sphere. 
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The technical and experimental advances associated to the investigation of micrometric 
and nanometric devices has boosted the research on dispersive forces, like Casimir and van 
der Waals forces which arise from fluctuations of the quantum vacuum due to bound- 
aries or external fields. As it is well known, Casimir |2| predicted an attractive force between 
two perfect conductors parallel plates separated by a given distance. The Casimir theory 
was extended latter to dielectric materials by Lifshitz in 1956 jjj. Due to the experimental 
difficulties inherent in handling two parallel plates at small distances, Derjaguin and col- 
laborators measured the Casimir force between a glass plate and a relatively large spherical 
lens 0| and they developed the so called Dejarguin approximation or proximity theorem. 
This theorem allows the estimation of the Casimir force per unit area between a sphere 
and a plate if the energy per unit area between two parallel plates is known. However, this 
theorem only holds when the radius R of the sphere and the distance z between it and the 
substrate satisfies the relation, z/R -C 1. Most of the accurate measurements of the Casimir 
forces perked in recent times re ly upon this t h eore m flflQ. However, it is not dear up 
to what limit this approach remains valid, and of course, for R ~ z, an exact calculation of 
the Casimir force between a sphere and a planar surface becomes essential. 

Within the Casimir effect context, the problem of a particle over a substrate reduces to 
finding the effects of the environment acting on the particle. This is not a simple problem 
since it must be solve in a self-consistent fashion, since the environment, in this case the 
substrate, has an effect on the particle and the particle also affects the environment. In a 
previous work the optical properties of a spherical particle lying above a substrate have 
been determined through its response to a local electromagnetic field. In this case, this field 
is produced by an external or applied field, plus a field coming from the presence of the 
substrate. The latter was introduced in the problem as a charge distribution induced on 
the particle due to the presence of the substrate which modifies the response of the particle 
itself. This procedure is implemented by considering the interaction of the induced dipole 
moment of the particle with its own image dipole in the substrate. However, this approach 
fails to study the optical properties of nanoparticles that are very close to the substrate, 
such that the "image-dipole" model must be extended in order to include higher-order 
multipolar moments. Previous extensions of the "image-dipole" model have shown that the 
importance of these multipolar excitations increases with the proximity of the particle to 
the substrate Q]. This "image-dipole" model was also considered by Casimir and Polder, 
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in 1948, to find the London or van der Waals forces of a polarizable atom near a perfect 
conductor plane Q. They found a correction to the van der Waafs force by considering the 
influence of retardation. In 1998, Ford |10| also considered the "image- dip ole" model to 
calculate the Casimir Force between a perfectly conducting wall and a large sphere with a 
Drude dielectric function. In early experiments to measure the Casimir-Polder interaction, 
Shih and Parsegian [ljj measured the deflection of atomic beams by gold surfaces obtaining 
results that are inconsistent with the "image-dipole" model. The measurements of Shih and 



121, indicate that the Casimir 



Parsegian h|, as well as other experimental observations 
or van der Waal force between a polarizable particle and a planar substrate should involve 
more complicated interactions than the simple "image-dipole" model. 

In this paper, we calculate exactly the dispersive force, or Casimir force in the non- 
retarded limit, between a nanosphere and a plane. The force is obtained considering all the 
high-multipolar excitations in system, that are shown to enhance the force by several orders 
of magnitude as a nanosphere approaches the substrate. We use a spectral representation 
(SR) method jl^ to calculate the force, that yields the frequencies (u s ) of the proper electro- 
magnetic modes of the sphere-substrate system trough a diagonalization process. Therefore, 
the zero-point energy is easily calculated as 8 = J2 S V^^V Another advantage of the SR is 
that it separates the contribution of the dielectric properties of the sphere from the contri- 
bution of its geometrical properties. Both properties of the SR allows to perform systematic 
studies of dispersive forces in this system. 

We consider a nanometric-size sphere of radius R located at a a minimum distance z(< R) 
from a substrate. In this case, the radius of the sphere and the minimum separation between 
the sphere and the substrate, are smaller than the characteristic length of the system |l^ |. 
therefore, it is not necessary to consider retardation effects. The fluctuations of the electro- 
magnetic field induce a multipolar moment which also induces a charge distribution in the 
substrate that affects the multipolar moment on the sphere, trough a self-induced polariza- 
tion process. Then, the lm-th multipolar moment on the sphere, Qi m , is induced by a local 
electromagnetic field, V^ c , which is the sum of the induced field due to the quantum vacuum 
fluctuations at the zero-point energy (V^™ c ), plus the induced field due to the presence of 
the substrate (V^ b ). Therefore, the Qi m of a sphere with polarizabilities a.i m {uj) is given 
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by [1 



Qlm = - —J- ^ Oil m (uj)V^ C , (1) 
(21 + 1) . . r T rvac i T/subl 

= ^ — aim(uj) [V lm +V lm \, 

where the Im-ih multipole moment induced in the sphere is defined as 

Qim= I r n p(r')Y lm (0',<p')d 3 r', (2) 

where p(r') is the charge density in the sphere, Y\ m are the spherical harmonics, and the 
integral is performed over the volume v s of the sphere. The /m-th component of the local 
field can be written using a multipolar expansion [15] that yields to 
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where Qv m > is the I'm' -th induced multipolar moment in the substrate which is located at 
r = (2(z + R), 6 = 7r, ip) from the center of the sphere, and A l l™' is the matrix that couples 
the interaction between the multipolar distribution on the sphere and substrate The 
induced Z'm'-th multipolar moment in the substrate or the "image-multipole" moment is 
related with the Qi m on the sphere by 

Ql'm' = {-l) l ' +m 'fc{u)Ql> m >, (4) 

where f c (uS) is a contrast function that together with the factor (— l) z '+ m ' allows to satisf 
the boundary conditions of the electromagnetic field on the sphere-substrate system 
Here, we have assumed that the substrate is electrically neutral and has a local dielectric 
function e sub (u;), and the ambient between the sphere and substrate is the vacuum. Let us 
suppose that the sphere is homogeneous, and has a local dielectric function e sp h(o;), such 
that it has azimuthal symmetry, and its polarizabilities are independent of m, and are given 
by H 

«iH = n^T^ R 2l+1 - (5) 
lK ' l[e sph (u) + 1] + 1 K J 

To find the proper electromagnetic modes of the sphere-substrate system, we need to solve 
Eq. (3) by a self-consistent procedure or by inverting L(2L + 1) complex matrices of L x 
L dimension, where L is the largest order of the multipolar expansion. To avoid these 
cumbersome procedures, we use instead a spectral representation that allows to find the 



proper electromagnetic modes trough the diagonalization of a real and symmetric matrix. 
The discussion of the above method can be found in Ref. 13. 
First, we rewrite the sphere polarizabilities from Eq. (JHJ) as 

= ~r\-R 2l+ \ (6) 
n — u{uj) 

where u(uj) = [1 — e sp h(^)] _1 and n! = 1/(21 + 1). One should notice that the poles of 
Eq. (JHJ), that is u{uji) = n l , yield the frequencies of the proper electromagnetic modes of the 
isolated sphere, therfore u(uj) is known as a spectral variable Q]. Now, we rewrite Eq. (jHJ) 
as an eigenvalue equation as 

[-«(">) V + H t] V = 9» (?) 

where we have simplified the notation by writing fi = (/, m), and 

_ Qim _ (IR 21+1 ) 1/2 ac , 

M (lR 2l+i yi 2 ' M 47r m ' 

and 

H( = n* V + / c ( W )^^(/0 1/2 <- (9) 

It is shown in Ref. 15 that Afi is a symmetric matrix that depends on the distance between 
the "image-multipole" and the sphere as, l/[2(z + R)] l+l>+1 , such that H£ is dimension- 
less and symmetric, and only depends on the geometry of the system through the ratio 
z/R. Consider the case when the contrast factor function f c (u) is real 
also real, and we can always find a unitary transformation that diagonalizes it, such that 
Eiiu'(^)" 1 ^^ = ^ n u^uu'i being n v the eigenvalues of the matrix H£ . Then, the 
solution of Eq. (JZj) is given by 

^ = -E G mV, (10) 
**' 

where is a Green's operator, whose \x\i! element can be written in terms of the Ujf and 
the eigenvalues of Eq. © , as 

jjml (Tjm"l" \-l 
nt 1 f„A — nml _ u m"l"\ u m'V I 

<VW - ^m'V ~ 2^ " _ „l" • ^ li > 

m"l" m " 

The poles of G^, lf (u), that is u(u l m ) = n l m , yield the frequencies of the proper electromagnetic 
modes of the system. If the sphere is isolated there is no substrate, therefore, the dielectric 
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then Htt is 



function of the substrate is replaced by the one of the vacuum and f c (uj) = 0, such that 
n l m = raj,. 

We now calculate the interaction energy as the difference £{z/R) = £ Z / R — Soo, where 
£ z /r is the energy when the sphere of radii R is at a distance z from the substrate, and £oo 
is the energy when the sphere is far from the substrate (z — > oo) or isolated, that is, 

£(,/fl)=£^-£M. (12) 

lm I 

The results presented here are calculated as follows. First, we construct the matrix H for 
a given z/R, choosing a maximum value of multipolar excitations that ensures convergence 
in the eigenvalues and eigenvectors of H, and then, we diagonalize it. We consider an 
explicit expression for the dielectric function of the sphere, and then we calculate the proper 
electromagnetic modes trough the relation u(u l m ) = n l m . Once we have u l m , we calculate the 
energy according with Eq. (jl2J) . The largest order of the multipolar expansion considered 
in this work was L = 2000 at the minimum distance of z/R = 0.001, because multipolar 
modes of higher order than L do not contribute to the energy S(z/R) at such distance. 
Here, we use the Drude model for the dielectric function of the sphere, therefore, u(u) = 
[u(u + i/r)]/ujp, where u p is the plasma frequency and r is the relaxation time. We present 
results for aluminum (Al) spheres with ftw p = 15.80 eV, and (rtUp) -1 = 0.04. The substrate is 
sapphire whose dielectric function is real and constant in a wide range of the electromagnetic 
spectrum, with e su b = 3.13 and contrast factor f c {uj) = —0.516. 

In Fig. 1, we show the energy calculated according with Eq. ()12|) as a function of z/R. 
We show results when all multipolar effects are taken into account (solid line), as well as 
results when only quadrupolar effects (dashed line), and dipolar effects (doted line) are 
considered. In general, we observe that the energy shows a power law of (z/R)~ 3 when only 
dipolar effects are considered while for the curve corresponding to quadrupolar multipoles 
shows three different regions: (i) when z/R > 7 the energy shows a power law of (z/R)~ 3 , 
that means that only dipole-dipole interactions are important, (ii) when 2 < z/R < 7 the 
energy shows a power law of (z/R)~ 4: , indicating that dipole-quadrupole interactions are 
important, and (iii) when z/R < 2 the energy shows a power law of (z/R)~ 5 , such that only 
quadrupole-quadrupole interactions become important. 

At small distances, the energy of the system can increase up to three orders of magnitude 
when all the contributions from high-multipolar moments are considered, as compared with 
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the energy when only dipolar or quadrupolar moments are taken into account. Notice 
that when we are calculating the energy with all the high-multipolar contributions, we are 
solving the problem exactly Then, the exact energy curve shows a power law of (z/R)~ 3 
when z/R > 7, when 2 < z/R < 7 the energy shows the same behavior as the curve where 
quadrupolar effects are included, and at smaller distances (z/R < 1) the energy increases 
sharply as the minimum distance also does, and it is not possible to assign a power law. This 
means that high-multipolar contributions become very important as the sphere approaches 
the substrate. In other words, the induced field in the substrate becomes very inhomogenous 
when the sphere is almost touching the substrate and then the attractive force between the 
sphere and the substrate increases tremendously. 

In Fig. 2, we show the Casimir force for an Al sphere over sapphire calculted as F = 
—dS/dz. We show the force when all multipolar interactions are taken into account, as well 
as up to dipolar, and up to quadrupolar interactions are considered. In agreement with 
the results for the energy, we observed for the force that multipolar effects become evident 
when the minimum distance between the sphere and the substrate is smaller than R. As the 
sphere approaches the substrate, the attractive force suddenly increases up to four orders 
of magnitude as compare with the dipolar interactions. At large distances (z > 2R) the 
force can be obtained exactly if only up to quadrupolar interactions are considered. We 
also obtained that for z > 7R the interaction between the sphere and the substrate can 
be modeled using the dipolar approximation [~L3|. In Fig. 3, we show the Casimir force for 
spheres of different radii, where all the high-multipolar modes were taken into account. We 
observed that when the sphere is almost touching the substrate, the attractive force is larger 
for the sphere with smaller radii, however, the force decays faster at larger distances when 
the sphere becomes smaller. Then, we observe that at a minimum distance of z = 50 nm, 
the force for the sphere with R = 100 nm is ten times larger that the force for the sphere 
with R = 10 nm. Notice that at such distance only dipolar interactions are important if 
R = 10 nm, while multipolar contributions are of importance if R = 100 nm. 

In conclusion, we have shown that the dispersive or Casimir force is enhanced by several 
orders of magnitude as a nanoparticle approaches a substrate. We have calculated exactly the 
interaction energy between the nanoparticle and the substrate using a spectral representation 
formalism. We found that as the separation between the particle and the substrate decreases, 
more multipolar moments contribute to the energy. On the other hand, for large separations 
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the dipolar term is the dominant one, like in the Casimir and Polder model. The radius of 
the nanoparticle plays an important role, since at small distances the force increases if the 
radius decreases, while at large distances the contrary occurs. The increment of the force at 
small separations could explain the physical origin of the large deviations observed in the 
deflection of atomic beams by metallic surfaces. However, specific experiments have to be 
perform to prove the this. 

This work has been partly financed by CONACyT grant No. 36651-E and by DGAPA- 
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FIG. 1: Casimir energy as a function of the ratio z/R. 
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FIG. 2: Casimir force as a function of the minimum distance z for a sphere with R = 50 nm. 
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FIG. 3: Casimir force as a function of the minimum distance z for spheres with R = 10, R = 50, 
and R = 100 nm. 
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